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Elementary particles possess quantized values of charge and internal angular momentum or 
spin. These characteristics do not change when the particles interact with other particles or fields 
as long as they preserve their entities. Quantum theory does not explain this quantization. It 
is introduced into the theory a priori. An interacting particle is an open system and thus does 
not obey conservation laws. However, an open system may create dynamically stable states with 
unchanged dynamical variables via self-organization. In self-organized systems stability is achieved 
through the interplay of nonlinearity and dissipation. Can self-organization be responsible for 
particle formation? In this paper we develop and analyze a particle model based on qualitative 
dynamics and the Feigenbaum universality. This model demonstrates that elementary particles can 
be described as self-organized dynamical systems belonging to a wide class of systems characterized 
by a hierarchy of period-doubling bifurcations. This semi-qualitative heuristic model gives possible 
explanations for charge and action quantization, and the origination and interrelation between 
the strong, weak, and electromagnetic forces, as well as SU(2) symmetry. It also provides a basis 
for particle taxonomy endorsed by the Standard Model. The key result is the discovery that the 
Planck constant is intimately related to elementary charge. 



I. INTRODUCTION 

Our world appears to be quantized. The most basic 
creatures of matter/fields, the elementary particles, pos- 
sess fixed values of strong, weak, and electromagnetic 
charges and angular momenta. These fixed dynamical 
variables do not change when a particle interacts with 
the surrounding vacuum and with other particles/fields. 
In classical mechanics, two different classes of systems 
can exhibit stability: the conservative systems and the 
dissipative systems [2J [13] . 

In conservative systems, dynamical variables are con- 
served due to the existence of symmetries (the Noether 
theorem) . The finite motion of a classical conserved sys- 
tem can be described by a closed loop trajectory that 
is parameterized by the corresponding conserved vari- 
able. Different trajectories densely fill the state space, 
and an infinitely small perturbation can shift the system 
from one trajectory to another. The new trajectory cor- 
responds to a different value of the conserved variable 
and is as "stable" as the previous one. This system is 
not asymptotically (absolutely) stable in the sense that 
a small perturbation does not asymptotically fade out 
after the interaction. Thus, conservation is conditional 
and requires the conservative system to be closed (i.e. 
it cannot interact with the rest of the world). However, 
elementary particles interact with external fields; they 
are open and the stability of their dynamical variables 
cannot be explained by the conservation laws of classical 
mechanics. 

In contrast, asymptotic stability is quite common in 
nonlinear dissipative systems and manifests itself as a 
phenomenon of self-organization [TJ]. The theory of 
dynamical systems, which includes nonlinear and dissi- 
pative phenomena, has succeeded in understanding the 
origin of numerous patterns like vortices, domain walls, 



pinches, various sorts of waves, Benard cells, linear and 
point defects, etc. that occur in dynamical media such as 
fluids, the atmosphere, chemical reactions, gaseous and 
solid-state plasmas, laser cavities, electric circuits, cellu- 
lar automata, etc. 

Self- organized systems (SOSs) are open. Their stabil- 
ity comes from the interplay between nonlinearity and 
dissipation. The stability is dynamical and is achieved in 
states that are far from thermodynamic equilibrium. In 
many cases, SOS phase portraits represent trajectories 
spiraling toward closed loops called attractors. Attrac- 
tor closure implies stability in the corresponding dynam- 
ical variables. Near the attractor an SOS behaves like a 
conservative system. However, the SOS stability is not 
conditional. Due to dissipation it is asymptotic — small 
perturbations fade out with time and the system returns 
to its original attractor [2j 1X3] . 

Self-organization is a ubiquitous phenomenon and has 
been observed at different scales of matter. The observ- 
able part of the Universe looks like a hierarchy of self- 
organized structures, starting from galaxy super-clusters 
to stars and planets to atmospheres and ecosystems to 
living organisms and their organs to cells and microor- 
ganisms, and all the way down to molecules and atoms. 
All of these systems are far from thermodynamic equilib- 
rium they are dynamical systems. Hubble's law demon- 
strates the dynamical state of the visible universe at the 
largest scales. Stars are born and die. Novas explode. At- 
mospheres seethe. Plants grow. Animals breathe. Cells 
self-reproduce. Even vacuum exists in a state of ther- 
modynamic non-equilibrium — it is filled with CMB radi- 
ation, neutrinos, and other excited fields. We have all 
reason to assume that this vibrant multi-level dynami- 
cal pyramid interlaced with self-organization can be ex- 
tended to the subatomic level. But, can the phenomenon 
of self-organization explain the stability of elementary 
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particles and the quantum nature of their internal angu- 
lar momenta, charges, and masses? 

Traditionally, particle physics belongs to the frame- 
work of quantum theory (QT). QT acknowledges parti- 
cle openness and the fact that particles interact with the 
surrounding medium. QT utilizes Lagrange-Hamiltonian 
mechanics of conservative systems but in a different way 
than its classical counterpart. The major difference here 
is the introduction of special constraints, the quantiza- 
tion rules. Quantization provides absolute stability to 
the particles, a stability that is absent in classical me- 
chanics of conservative systems. In QT, dynamical char- 
acteristics such as charge and action are postulated to 
have only discrete or fixed values, and these constraints 
dictate which dynamical trajectories are permitted and 
which are forbidden. The permitted trajectories consti- 
tute a discrete set in the corresponding state space. Now 
a small perturbation cannot shift a particle from its per- 
mitted trajectory to a nearby trajectory for arbitrarily 
long periods of time because the latter trajectory is for- 
bidden. Thus, the perturbation is expected to die within 
the time interval limited by the Hciscnbcrg uncertainty 
principle. Such processes, known as virtual processes, 
strongly resemble dissipation, and the most probable QT 
"trajectories" , the eigenstates, strongly resemble SOS at- 
tractors. The possible importance of dissipation in the 
foundations of QT has been discussed by 't Hooft [15 . 

Despite being based on the superposition principle, 
QTs framework includes elements of nonlinearity. The 
commonly used perturbation theory is concerned with 
nonlinear corrections. The renormalization technique 
assumes that charges (coupling constants) and masses 
depend on the perturbation level. Non-Abelian Yang- 
Mills theories of weak and strong interactions use charged 
bosons that interact with one another and make the non- 
linearity even more profound. Abrupt transitions from 
one state to another as well the entire measurement pro- 
cess (the so-called collapse of the wave function) are 
"super-nonlinear'' they arc discontinuous. The latter 
also breaks time symmetry, and is thus a dissipative pro- 
cess. 

We can see that QT objects possess some important at- 
tributes of self-organized systems: nonlinearity and dissi- 
pation. The QED triumph (one example being the strik- 
ing accuracy in calculating the electron magnetic mo- 
ment) would be impossible if an electron were treated 
strictly as a closed linear conservative system. Despite 
the similarity between QT objects and SOSs, absolute 
stability emerges in different ways in the two frameworks. 
In SOSs it is an outcome of the theory. In QT it is intro- 
duced a priori via quantization rules and fundamental 
quantum constants like h, e, a, sm(9\y), sin (9c), etc. 

In this paper we build a model describing elemen- 
tary particles as self-organized systems. Our approach is 
heuristic and phenomenological. We are not concerned 
with a specific form of differential equations. Dissipa- 
tive systems are immune to small perturbations and their 
"global" behavior is not very sensitive to the details of the 



governing equations. To understand the overall SOS be- 
havior, it is often sufficient to use qualitative analysis and 
study a prototype-system that belongs to the "proper" 
dynamical class. 

In the Standard Model (SM) we find a number of differ- 
ent doublets such as particle anti-particle, spin-up-spin- 
down, proton-neutron, lepton-quark, u-quark-d-quark, 
electron-neutrino, etc. We can even organize all the fun- 
damental fermions as a doublet pyramid resembling a 
"phylogenetic" tree, where each fork represents a new 
doublet (FigjTJ). In many cases these doublets are seen 
as two states of a single particle. The transformation 
from one state to the other can be viewed as a rotation in 
some complex internal space. A remarkable property of 
these rotations is their SU(2) symmetry, which requires 
rotation through 4-7T rather than the usual 2n to return 
to the original state. This period doubling and the bifur- 
cating structure of the "phylogenetic" tree astonishingly 
resemble the period-doubling bifurcation diagrams that 
are often found in systems with non- linear dynamics, and 
are typical for a large class of dissipative systems. The 
similarity is clearly evident from the comparison of Fig[l] 
with Fig(2] which shows the period-3 stability window 
from the logistic map bifurcation diagram. The three 
branches in the bifurcation diagram resemble three fam- 
ilies in the particle diagram. 

Our observation suggests that particles may belong to 
a class of systems in which the dynamics are character- 
ized by period-doubling bifurcation diagrams. The class 
is large, and includes a number of dynamical systems 
that not only possess similar bifurcation diagrams, but 
even share the same scaling properties. This universality, 
originally discovered by Feigenbaum [3 [6j |7] , allows us to 
add to our qualitative analysis some quantitative results 
that can be obtained when we represent an elementary 
particle with the universal Feigenbaum function or with 
another function from the above-mentioned class. 



II. BUILDING A MODEL 

Let us take the electron as an example. To visualize 
the complexity of the internal electron dynamics, let us 
imagine a hypothetical situation that can arise as a re- 
sult of a negatively charged fluctuation in vacuum. The 
fluctuation polarizes the surrounding vacuum and cre- 
ates a positively charged halo around itself. The halo 
lowers the local electric potential, and the original fluc- 
tuation becomes denser and more confined. This in turn 
affects the halo. The positive feedback described above 
competes with the negative feedback caused by charge 
diffusion and self-repellence. Thus, the symmetric halo 
becomes unstable and breaks up into separate positively 
charged fragments, which then create secondary nega- 
tively charged halos around themselves. This process re- 
peats itself ad infinitum at smaller and smaller scales. 
Moving fragments create currents and magnetic fields, 
adding more complexity to this turbulent system. Mov- 
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FIG. 1: Three fermion families can be organized as bifurca- 
tion diagrams. 



FIG. 2: Period-3 window from a bifurcation diagram of the 
logistic map. 



ing positive and negative fragments experience attractive 
and repellant forces from their neighbors. Both are pro- 
portional to the fragment charges and currents. Due to 
the finite distance between fragments and the limit of 
the speed of light, feedback is delayed. Depending on the 
strength of the perturbation, the original fluctuation ei- 
ther relaxes to a state of thermodynamic equilibrium or 
the system bifurcates to another state where the delayed 
interplay between attraction and repulsion gives birth to 
a dynamically stable spatiotemporal pattern. We suggest 
that an electron is one such self-organized system. 

Due to the fractal structure and openness of this sys- 
tem, there can be infinitely many interactions among the 
particle parts. To describe these dynamics we need an 
infinite-dimensional state space. The complexity of this 
situation resembles a gas ensemble, or rather, a turbu- 
lent fluid. A practical approach to this type of complex 
system is to find a few collective variables such as tem- 
perature, concentration, pressure, or convection velocity, 
that effectively describe the dynamical state of the sys- 
tem in a low-dimensional state space |13j . In our case, 
it is most natural to assign this role to charge. Thus, 
we consider a low-dimensional state space as a nonlin- 
ear vector field where q is the effective collective 

dynamical variable, and attempt to understand why q is 
quantized. In QT, electrical charge is the coupling con- 
stant between the particle and the external electromag- 
netic field. Since our model is primarily concerned with 
internal dynamics, we also assume that charge is the cou- 
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pling constant — albeit a running coupling constant — in 
the interaction between the particle and its internal field, 
the field that is created by the particle itself. In other 
words, charge defines the feedback that is responsible for 
the particle's self-organization. 

We assume that as an SOS, an electron possesses 
asymptotic stability. It exercises finite motion in the 
corresponding state space, which can be represented by 
dynamical trajectories in that state space. Dissipation 
ensures that such trajectories spiral toward their limit 
cycles, the attractors. We assume that the attractors 
describing dynamical equilibrium can be parameterized 
by charge q. To make this parameterization sensible, we 
also assume that the system possesses some "inertia" , 
which means that external perturbations, though capa- 
ble of changing the q- value (q is a running parameter), 
cannot do it abruptly. Thus, q is preserved for at least 
a few cycles, obeying the so-called adiabatic constraint 
[S]. If the dynamics are not chaotic, the attractors are 
closed curves as in the case of conservative systems. We 
can use this similarity to describe the dynamics in terms 
of generalized action-angle variables [5] , where angle ip is 
a cyclic or ignorable coordinate, and the generalized mo- 
mentum or reduced action (or simply action), J = S/2ir 
(S is the total action accumulated during the entire cy- 
cle), is a "conserved" variable. The latter can also be 
used to parameterize the attractors. Having selected q 
as a dynamical parameter, we would like to connect the 
J-parameterization with the ^-parameterization. This 
"canonical" transformation can be accomplished via di- 
mensional analysis if we define J as: 

J = m\ (i) 

where 7/ is a conversion constant that has the physi- 
cal dimensions of electromagnetic impedance. To be 
more specific, we assign r\ to be the vacuum impedance, 
V = \J Mo/eo; where eo is vacuum permittivity and fiQ is 
vacuum permeability. 

To further simplify our analysis, we replace each con- 
tinuous trajectory in the state space with a set of points 
ipi (one point per loop) that are selected using a proce- 
dure called the Poincare section [T2] (Figj3J. 

After applying the Poincare section, we obtain a one- 
dimensional recurrent map called a Poincare map: 

V7 fe = F ($k-i) , (2) 

where F Otpj is a recursive function. According to our as- 
sumptions, this function can be parameterized by charge 
q and denoted as F q yip^J , or by action J and denoted as 

Fj (V). 

Based on similarities between the particle "phyloge- 
netic" tree diagram (FigJIJ and the period-doubling bi- 
furcation diagram (FigpJ, we assume that the dynamical 
electron is similar to systems in which the dynamics pos- 
sess a period-doubling bifurcation structure. The class 




FIG. 3: Selection of map points ipi using the Poincare section. 



of period-doubling bifurcation systems is wide [4 . It in- 
cludes all maps with smooth unimodal (having a single 
extremum) recursive functions (see examples in Figj4a-c), 
and many known SOSs such as mechanical and electronic 
oscillators, Benard cells, Belousov-Zhabotinsky chemical 
reactions, Couette-Taylor flow, etc. 

A bifurcation diagram maps limit points Xoo (vertical 
axis) against control parameter values (horizontal axis). 
It represents a hierarchy of period-doubling bifurcations. 
Each bifurcation is a phase transition that is accompa- 
nied by a doubling in the number of limit points (at- 
tractor loops). The distance between bifurcations pro- 
gressively shrinks with parameter, and the scaling fac- 
tor quickly converges to a number called the Feigenbaum 
constant. Our three example maps (Figj4^,-c) have dif- 
ferent recursive functions, yet their bifurcation diagrams 
still look very similar. They are called quadratic maps 
because the extrema of their recursive functions can be 
approximated by quadratic parabolas. According to the 
Feigenbaum universality [5J O |7], the scaling factor for 
these maps, as for all unimodal quadratic maps, is in- 
dependent of all other details in their recursive func- 
tions, and always converges to the Feigenbaum num- 
ber 5 — 4.669.... We assume that the "electron" map 
represented by equation Q also belongs to the class of 
quadratic maps. The corresponding hypothetical bifurca- 
tion diagram of the particle recursive function Fj is 

shown in Fig|4ji. We select J (not q) as the parameter for 
the following reason. In all quadratic maps, parameter A 
controls the strength of the feedback or self-interaction. 
In electromagnetic theory, feedback is a two-step pro- 
cess: 1) the excitation of the electromagnetic field by the 
charge distribution, and 2) the formation of the charge 
distribution by this field. The "intensity" of each step is 
proportional to q, and the "intensity" of the entire feed- 
back is proportional to q 2 , or assuming equation |lj, to 
J. Thus, J plays the same feedback role in electron self- 
interaction as parameter A does in quadratic maps. For 
convenience, we place the origin of the J-axis, J = 0, at 
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FIG. 4: Period-doubling bifurcation diagrams for (a-c) three different quadratic maps and (d) a hypothetial diagram for 
particles. 



the Feigenbaum point [12], and the origin of the ^-axis, 
ip = 0, at the point where the electron recursive function 

reaches its hypothetical extrcmum. 

Bifurcation points divide maps into segments of sta- 
bility [IS]. Each map has a continuous interval where 
it converges to a single fixed point that corresponds to 
a single loop attractor (we do not account for the de- 
generated case of thermodynamic equilibrium where the 
attractor is just a point). In the next parameter inter- 
val, the original map looses stability but its second it- 
eration ipk = Fa (^Fa (^Pk-ij^j still converges, implying 
the existence of a period-2 attractor. Now we have two 
limit points in the bifurcation diagram K^oo) ^io) > which 

means that the attractor crosses the Poincare section 
two times, and represents a closed double-loop trajectory 
(Fig(5|. Period-doubling bifurcation manifests a phase 
transition that changes the structure of attractors pQ. 
It can be considered as the emergence of a second peri- 
odic motion with a doubled period (schematically shown 
in Fig{5] by small circles) whose dynamical mode is or- 
thogonal to the original mode. Describing the new dy- 
namical state requires an additional degree of freedom. 
The description of a new trajectory requires two complex 
numbers, and taking into account period-doubling, can 
be represented by a two-component spinor. This can be 



considered as a doubling of the state space dimensionality 
by the bifurcation. 

As the control parameter reaches its second critical 
value, another period-doubling bifurcation occurs. The 
stability shifts from a period-2 attractor to a period-4 at- 
tractor. The state space experiences a similar metamor- 
phosis and its dimensionality doubles again. The projec- 
tion of the attractor onto the original state space now 
represents four loops and needs a four-component spinor 
for its description. At the next critical value a period-8 
attractor replaces the period-4 attractor and so on. This 
progression ends when the number of branches, the di- 
mensionality of the dynamical space, and the number of 
spinor components reaches infinity, the so-called Feigen- 
baum point. 

One quadratic map that has been extensively explored 
is now considered to be the paradigm for the entire class. 
It is called the logistic map (Figj4^i). We will use it in 
our numerical explorations. 

The logistic map is given by the equation: 

x k = Axk-i (1 - x k -i) , (3) 
where A is the control parameter. 



FIG. 5: Double-loop attractor can be represented by a two- 
component spinor. 



III. ORIGIN OF QUANTIZATION - AN 
INHERENT FEATURE OF SOS DYNAMICS 

Now let us examine how the dissipation rate varies 
with the control parameter in the stability intervals, i.e. 
between bifurcations. 

The dissipation rate D can be defined as the inverse 
number of steps k (iterations), D = 1/fc, required for 
the system to go from a randomly selected initial point 
xq to the vicinity of the corresponding fixed point, i.e. 
within the interval (xk=oo — £, Xk=oo + e) for any suffi- 
ciently small e. Fig|6] demonstrates the dissipation rate 
of a map given by equation Q for the first three stability 
intervals when xq = 0.367, koo = fciooch an d £ = 10~ 14 . 

We notice that Z)(A)-curves have a sharp maximum 
close to the center of the corresponding stability interval. 
These maxima occur when attractors have limit points 
that coincide with points where the recursive function 
reaches its extremum: x = 0.5 for the logistic map and 
ip = for the "electron" map. Here the stability parame- 
ter, called the Lyapunov exponent and determined as A = 
{I /n) hx{\dF (x ) /dx\ \dF ( Xl ) /dx\ ...\dF (x n ) /dx\) 
|12j . has a singularity, A — > — oo. The corresponding 
attractors are called supercycles or superattractors. 
Because superattractors play an important role in our 
model, we denote their parameter values with subscripts 
representing the stability interval number. For example, 
Ji and qi are the action and charge corresponding to the 
superattractor of the i-th stability interval. 

Let us conditionally divide each converging trajectory 
into two parts, the first being the transient spiral, and 
the second being the vicinity of the attractor. For each 
selected trajectory, the higher the convergence rate, the 
shorter the transient time, the more time the system 
spends near the attractor, and thus the higher the proba- 
bility of finding the system in the vicinity of the attractor. 

Despite the adiabatic constraints, J is still a free run- 
ning parameter. Vacuum fluctuations or other external 
noise can kick the electron from one state trajectory to 
another and from one J-value to another. However, due 
to the profound differences in dissipation rates for differ- 
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FIG. 6: Dissipation rate as a function of the control param- 
eter for three intervals of stability (bottom) and the corre- 
sponding segments of the bifurcation diagram (top). 

ent attractors (Fig(6]), the average time spent near the 
superattractor is much longer than the time spent near 
other attractors, i.e. the superattractors are by far the 
most "attractive" attractors. Hence, the probability of 
finding the system near a superattractor is much higher 
than the probability of finding it in other regions of the 
state space. Thus, we come to the important conclusion 
that in our self-organized system there exists a set of 
special, discrete, most probable, and super-stable closed 
trajectories. This can be interpreted as quantization, a 
phenomenon similar to the one postulated in QT. The 
difference is that we did not introduce this quantization 
a priori. It is an inherent feature of SOS dynamics. The 
quantization of attractors imposes a quantization of their 
parameters. In reference to our model, this means that 
both the action J and the charge q have preferred discrete 
values, i.e. they are also quantized. 

It is interesting that in order to find superattractors we 
simply need to look for trajectories where the fixed point 
■0OO corresponds to the extremum of the recursive func- 
tion F (^ooj = F extr emumi which can be understood as 
a sort of variational principle. This variational principle 
can be extended to the situation when a particle interacts 
with an external field and the interaction deforms the re- 
cursive function, shifting its extremum and the extrema 
of its iteratives to other fixed points ipoo ( "eigenstates" ) . 
The deformation may also change the values of J n and 
q n . In the case of stronger interactions, the system jumps 
to the next stability interval with a different state space 
that has a doubled dimensionality and a different topol- 
ogy- 



IV. NUMERICAL SURPRISES 

According to our model an electron is a self-organized 
system for which the internal dynamics can be repre- 
sented by a discrete set of orbits in a state space or by 



7 



levels qi and Ji in the charge/ action scale. The specific 
values of charge and action depend on the level of excita- 
tion. The Feigenbaum universality makes it possible to 
obtain some quantitative relations between these param- 
eters at different excitation levels. The parameter values 
Jfe obey the same scaling law as the entire bifurcation 
tree — they have the same asymptotic behavior and con- 
verge to the same Feigenbaum delta [5J [7]: 



lim (Jk-i/Jk) = 5. 



(4) 



From equations Q and Q , the ratio between two charge 
values at two adjacent levels converges as: 



lim (g fc _i/g fc ) = 5 1 ' 2 . 

k—>oo 



(5) 



The convergence rate is usually high, and even for rela- 
tively small numbers k we can substitute equations Q 
and ([5]) with their approximations: 



•h-i/Jk = s 



and 



(Ik 



:-l/qk = ^ 2 . 



(G) 



(7) 



There is a special qi value that we refer to as q e . It cor- 
responds to the experimentally measured electron charge 
e: 



e= 1.6021... x 10~ 19 C*. 



(8) 



Let us explore the partice dynamics at different excita- 
tion levels around q e . As part of the exploration we will 
compare the dynamical variables Ji and qi at these levels. 

First we compare the parameters at the level where 
q e = e with the parameters two levels above where q e -2 = 
Se. Using equations ([I}, ([6|, and ^ we find that: 

Je-2 = 5 2 J e = r]S 2 e 2 = 2.1083 . . . x 10~ 34 J • s, (9) 

which equals twice the value of the Planck constant, 
2h = 2 (1.05457 ... x 10~ 34 J • s) , accurate to 0.04%. 
This implies that there is a direct connection between 
the two fundamental quanta e and h: 



= vVo/eo ( Se f / 2 - 



(10) 



Essentially ignored by QT, this relation plays a crucial 
role in our proposed model. The relation described by 
equation ( 10 ) carries important implications: 



1 . It indicates that quantization of charge and quan- 
tization of action have the same origin. 

2. The presence of the Feigenbaum delta implies the 
relevance of period-doubling bifurcation dynamics. 

3. The value of delta, 6 — 4.669..., suggests the 
involvement of dissipative dynamics (the period- 
doubling transition to chaos has also been found in 
Hamiltonian systems but in this case Sh ~ 8.721 . . . 



And finally, despite the fact that S is a fundamen- 
tal mathematical constant, it has only been used 
in one physical context, that of dynamical systems 
which experience transitions from smooth dynam- 
ics to turbulent dynamics, i.e. S belongs exclu- 
sively to chaos theory. The latter is essentially 
a non-quantum theory. For several fundamen- 
tal constraints (like superposition principle, uncer- 
tainty principle, and quantization itself), QT is in- 
capable of describing truly-chaotic systems. (The 
confusing term "quantum chaos" relates not to real 
chaotic systems, but rather to non-chaotic quantum 
systems whose classical counterparts are chaotic 
[3j [12].) Therefore, it is unlikely that one will be 
able to understand equation (10) (and the roots of 



charge and action quantization) within the frame- 
work of quantum theory. 



Equation ( 10 1 leads us to several other striking results. 
The fine structure constant, a — e 2 /Airsoch, can be ex- 
pressed exclusively via the mathematical constants 7r and 
6: 



= 2ir5 2 



1 

137' 



(11) 



By examining the next level, where charge is q e -\ = 
(5 1 / 2 e, we find that the "fine structure constant" is a e -i = 
(27r<5) -1 « 1/29 [TU], and corresponds to the experimen- 
tal value of the weak coupling constant. 

One more level above, where charge q e -2 = Se, we find 
that the "fine structure constant" is a e ~2 = (27r) = 
0.16, and corresponds to the experimental values of the 
low-energy (~15 GeV) strong coupling constant a s ob- 
tained through JADE data [151 ITS]. 

Our findings suggest that electromagnetic, weak, and 
strong interactions obey a hierarchy that can be ex- 
plained by period-doubling bifuration dynamics as shown 
in Fig{7] The (^-relations between coupling constants (but 
not their explicit values) were previously observed by 
Goldfain 0. 

Here are some additional examples suggesting that the 
Feigenbaum delta may have a profound relevance to par- 
ticle physics. 

According to the Standard Model, the photon and Z- 
boson are described as mixtures of the electromagnetic 
B-boson and the weak W 3 -boson. The ratios of their 
relative inputs are characterized by the weak mixing an- 
gle 6w Experiments show that sin 2 (9 W ) « 0.23 [TD], a 
value close to S~ x (see also [5]). 

Another example is the mixing of quarks from different 
families. Here, the Cabibbo angle 9c plays the role of the 
mixing angle. According to experiment, sin (0c) ~ 0.23, 
which is close to S^ 1 (see also [8]). 

Finally, a proton and a neutron can be considered as 
two states of the same particle, the nucleon. One pa- 
rameter that characterizes the two particles is a dimen- 
sionless physical constant called the (7-factor. Experi- 
mental values of the g-factor for a proton and a neutron 
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FIG. 7: Transformation of the strong force into the weak 
force into the electromagnetic force via bifurcations. 



are +5.585 . . . and —3.826 . . . respectively, which are two 
"asymmetric" numbers. We can restore the "symmetry" 
if we represent them as g P:Tl rs 1 ± S, which is consistent 
to an accuracy of a few percent . 



V. SPIN - ANOTHER POSSIBLE 
CONSEQUENCE OF THE SOS MODEL 

The physical dimension of the reduced action J is the 
same as the dimension of angular momentum, and at the 
level q e ~2 its value is 2h, suggesting that J may be re- 
lated to particle spin. Even more suggestive is the SU (2) 
symmetry, which is immanent of spin- 1/2 particles and 
plays a profound role in period-doubling bifurcation dy- 
namics. We propose the following construction in order 
to incorporate spin into our model. 

According to our model particle dynamics can be de- 
scribed by a 2™-component spinor, where n is the exci- 
tation level. To some approximation, these 2 n degrees of 
freedom can be viewed as independent. Similar to a gas 
system where the average energy is distributed among 
all degrees of freedom, we conjecture that in our parti- 
cle system the average action is shared by all degrees of 
freedom. We call this democracy the "equal action dis- 
tribution" rule. We also assume that particles can be 
located at different excitation levels and have different 
numbers of loops depending on the level they inhabit. 

For example, at the lowest excitation level, an elec- 
tron possesses a one-loop superattractor. This level cor- 
responds to small perturbations where we can neglect the 
probability of positron generation and spin flipping. At 
this level we assign the charge as q e = e and the action 
as J e — r/e 2 . At the highest level of electron excitation 
there are four quasi-independent electron states — spin- 
up electron, spin-down electron, spin-up positron, and 
spin-down positron. They can be represented by a four- 
component spinor (recall the Dirac spinor) or a four-loop 
attractor (Figj8^). These dynamics emerge two levels 
above q e — e, where q e -2 — Se and J e ~2 = S 2 rje 2 = 2h. 
Applying the "equal action distribution" rule we find that 
each degree of freedom in the four-loop superattractor 
has action 1/4 (2h) — h/2, which corresponds to the elec- 
tron spin. 

Photons come in two distinct polarizations and pos- 
sess two degrees of freedom that can be associated with 
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FIG. 8: Schematic attractor topology for (a) the electron, 
(b) the photon, and possibly (c) the graviton. (d) Electron 
branch of the "phylogenetic" tree. 



two-loop attractors (Figj8p). To find the lowest pho- 
ton level, we recall that a photon can be viewed as 
the fusion (annihilation) of two particles — the electron 
and the positron. In the particle "phylogenetic" tree, 
the electron and positron branches merge just one level 
above the lowest electron level (Fig{8]i). Here the charge 
is q e -i — 5 1 / 2 e. The two-loop photon superattrac- 
tor emerges at the next level where q e -2 — 5e and 
J e _2 = rj5J e -i = rj5 (5 1 / 2 e) = 2h. At this level the pho- 
ton action per degree of freedom is H, and corresponds to 
the photon spin. 

Continuing with the speculations, we come to the one- 
loop superattractor (Fig{8j:) at the level q e -2 = Se. Here 
action is J e _2 = Tj (5e) 2 = 2h and corresponds to a spin-2 
particle. The only elementary particle (albeit a hypothet- 
ical one) that is known to have spin 2h is the graviton. 

The bifurcation tree branch describing these three 
types of particles is shown in Fig j8ji. We refer to it as 
the "electron branch" . 



VI. PARTICLE ZOO 

Smashing particles in accelerators or other strong per- 
turbations may excite the system to levels beyond q e ~2- 
Such excitations may result in the emergence of more 
branches in the bifurcation tree (Figj9j). Each new bi- 
furcation doubles the number of degrees of freedom. As 
the system relaxes to its original level it is faced with 
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on. Adding new branches to the first family of particles 
assumes the existence of more bosons than is currently 
known from experiment. This oddity may be explained 
by the fact that it is difficult to distinguish among bosons 
of different families. Whereas fermions can be distin- 
guished exclusively by their mass (e.g. electrons, muons, 
and tau-particles have different masses), bosons are ei- 
ther massless or too heavy to be generated by our current 
particle accelerators. 
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VII. RELATION TO THE STANDARD MODEL 
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FIG. 9: Family I branch of the "phylogenetic" tree. 



more choices for particle self-organization. For exam- 
ple, after being excited to level g e _3 (Fig j9j , the sys- 
tem can relax to a lower level either through the elec- 
tron branch described in the previous section, or through 
a new branch (the neutrino branch), transforming into 
particles from a different group. The set of all options 
includes one more spin-2 particle, two more spin-1 parti- 
cles, and four more spin-1/2 particles. Altogether we now 
have four spin-1 bosons that can be associated with elec- 
troweak vector bosons, and eight fermions that can be as- 
sociated with spin-up/spin-down electrons/positrons and 
neutrinos / anti-neutrinos . 

The number of particles doubles again after an exci- 
tation to the level g e — 4- Now we have four spin-2 par- 
ticles, eight spin-1 bosons, and sixteen fermions. The 
bosons can be associated with gluons, and the eight ex- 
tra fermions with spin- up/spin-down u-quarks/d-quarks 
and their anti-particles. This particle zoo has a lot of 
similarities with the Standard Model taxonomy. How- 
ever, the analogy is not complete. For instance, some of 
the gluons coincide with electroweak bosons and photons. 

To accommodate particles from the second and third 
fcrmion families of the Standard Model we need to add 
more branches to our bifurcation diagram. If particle dy- 
namics belong to the period-3 window, this will be the 
end of the story. Otherwise, future experiments will dis- 
cover new particles belonging to a fourth family and so 



In our efforts we strive to be as close to the Standard 
Model as possible. Therefore, it is not surprising that 
the two models have a lot in common: SU (2) symme- 
try, the same types of forces, similar values of their cou- 
pling constants, the same number of fermions, electro- 
weak bosons, and gluons, etc. However, despite these 
similarities the two models have principal differences. 

The major distinction is that they are built on funda- 
mentally different premises. The Standard Model views 
conservation laws and the corresponding symmetries (or 
rather broken symmetries), both local and global, as 
the basic principles for particle existence. Our proposed 
model is necessarily based on dissipation, which inher- 
ently implies time-arrow asymmetry (the violation of T- 
invariance) and the existence of directional flows. These 
flows may also be responsible for violations of other 
global symmetries (like P-invariance, C-invariance, and 
CP-invariance). In addition, local symmetries and gauge 
invariance do not play the same fundamental and mys- 
tic role as they do in the Standard Model. The SU (2) 
symmetry simply reflects the relations among solutions 
of a special class of nonlinear equations (the recursive 
function and its iterations) that describe the evolution 
of particle dynamics under perturbations. Our model 
is formulated in a space-time independent framework. 
Both dynamical variables, action and charge, are Lorentz 
invariant. Thus, space-time symmetries are irrelevant. 
Some other profound differences between the two models 
are listed in the table in FigfTTJ] 

VIII. CONCLUSIONS 

We proposed a new framework for studying elemen- 
tary particles. This framework seems promising for pen- 
etrating into the previously inaccessible territory of the 
internal structure and dynamics of elementary particles. 
Particle self-organization is described via feedback that 
is controlled by running coupling constants (charges and 
actions). We explored the simplest system that mimics 
quadratic maps. We found sensible answers to a number 
of questions that have never been explained by quantum 
physics. We demonstrated that our model is capable of 
explaining the quantization phenomenon and has the po- 
tential to provide a basis for particle taxonomy and the 
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Standard Model 


Number and type(s) of 
fields 


Single nonlinear dissipative 

field with variable 
dimensionality and topology 


Set of independent fields 
obeying the superposition 
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Interplay between nonlinearity 
and dissipation 
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FIG. 10: Key differences between the proposed model and 
the Standard Model. 

relations among strong, weak and electromagnetic forces 
accepted by the Standard Model. 

Our key result is the discovery of the deep connection 
between two fundamental constants of quantum physics, 
the Planck constant h and elementary charge e. The rela- 
tion h = y/ no/eo (i5e) 2 /2 is experimentally verifiable and 
should be considered as an experimental fact irrespective 
of the correctness of our proposed model. Ironically, the 
two most fundamental quantum constants, h and e, are 
linked through the Feigenbaum S, a constant that be- 
longs to the physics of deterministic chaos and is thus 
exclusively non-quantum. 

Our results are assonant with 't Hooft's proposal that 
the theory underlying quantum mechanics may be dis- 
sipative [TS]. They also suggest that quantum theory, 
albeit being both powerful and beautiful, may be just a 
quasi-linear approximation to a deeper theory describ- 
ing the non- linear world of elementary particles. As one 
of the founders of quantum theory, Werner Heisenberg 
once stated, ". . . it may be that. . . the actual treatment 
of nonlinear equations can be replaced by the study of 
infinite processes concerning systems of linear differen- 
tial equations with an arbitrary number of variables, and 
the solution of the nonlinear equation can be obtained by 
a limiting process from the solutions of linear equations. 
This situation resembles the other one. . . where by an in- 
finite process one can approach the nonlinear three-body 
problem in classical mechanics from the linear three-body 
problem of quantum mechanics." |llj 
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